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OOM GARY 


Jeing Holzer's method of frequency calculation, the natural 
frequencies for the first two modes of torsional vibration of the 
wing were determined for a representative conventional airplane 
(p24-C, in the custonary manner, the fuselage being oonsidered as 
@ rigid body. Next, using a method developed by WN. 0. Myklestad 
of the Guggenhein Aeronautics Laboratory at the California In- 
stitute of Technology, conbined with Jolzer's method, the natural 
frequencies for the sane two modes of vibration were again deter- 
mined, but with the fuselage this tine being considered as flex- 
ible. 

A comparison of results of the two methods indicates that in 
considering the fuselage as being flexible, a decrease in the 
natural frequency of torsional vibration may be expected. For 
the particular airplane selected, this decrease enouted to 6.6890 
for the first mode of vibration and to 39.17% for the second. 

The investigation reported in this paper was entirely theo- 
retical and was perforned during the 1945-1944 school year at the 
Guggenhein Aeronautics Laboratory at the california institute of 
Technology, Pasadena, California under the direction and super- 
vision of Dr. N. 0. Myklestad, research associate in aeronautics 


at the Institute. 








INTRODUCTION 


In the a of aodern aircraft for higher and higher speeds, 
the designers are beconing increasingly more interested in the 
problens of flutter and vibration. One of these problems is that 
of the torsional vibration of the wings, which is dependent upon 
a nunber of factors, such as (1) the mass distribution both span- 
wise and chordwise of the wings thenselves and of all wnits sun- 
ported either on thea or within them, (2) the torsional stiffness 
of the wings, (3) the torsional! nonent anolied to the wings at the 
shifting center of pressure by the air loads, (4) the coupling be- 
tween the wings in bending and the wings in torsion, (5) the tor- 
Sional monent applied at the root of the wings by a flexing fuse- 
lage, and (6) the effect of conpressibility as looal velocities 
over the wing approach the velocity of sound. It is believed to 
be ule practice in the aircreft industry generally to consider 
all but the last two of the factors enumerated above. 

This paper then has as its objective the quantitative deter- 
Si aah ice of the effect on a representative large airplane of the 
fifth factor ome naweil above, nanely, the effect of the torsional 
monent avviied at the root of the wings by a flexing fuselage on 


the natural frequencies of the wing torsional vibration. 








DEPINITIO Cr SYWBOLS 


Fuselage aass concentrated at any station n, lbs. seconds 
squared per inch. 

Nunber of any wing or fuselage station. 

Mass monent of inertia of wing in.inch lbs. seconds squared 
or bending monwent of inertia of fuselage in inches to the 
fourth power. 

Convenient reference value of Te moment of inertia for 
fuselage as a whole. 

Modulus of elasticity of fuselage bending material in lbs. - 
per square inch. 

Angular deflection of wing at any station n in radians. 
Frequency of vibration in radians per second. 

Panel length of wing or fuselage between stations n and 
n+ | in inches. 

Shear at fuselage station n in lbs. 

Bending monent at fuselage station n in inch lbs. 


: f of fuselage tail at elastic axis. 


" tr? tt tf WW " 


nose 
Slope of fuselage axis at any station n. 


ff 8 " tail at elastic axis. 


tt t i te tt f 


nose 
Deflection of fuselage at any station n in inches. 


" " . tail at elastic axis. 








y, ~ Deflection of fuselage nose at elastic axis. 


ae Change in slope frm n to “+i due to a wit force at n. 


Ya, r n tt n n "  wea\ = "" f monent at . 
dp, - 8 " deflection frm 1 to n+! due toa unit force at “1. 
diy: - " % tt tt Y ¢t Y\ + \ 2] gy tf monent at n. 


ph - Slove of fuselage axis at extrene end of tail or nose. 
Ve | f, = Coefficients appearing in equation for &,. 
ae k t" 5 n n " Ya 
M. = Total coupling moment introduced into wing by fuselage at the 
elastic axis. 
= t/a ~ Non-dinensional syndol for fuselage bending monent of 
inertia at any station n. 


b= Opa, - 98, = Ens te - Non-dimensional symbol for inerease in 


fuselage bending acnuent of inertia fron 
station n to station n+!. 


& - Variable distance fron station n to any point in panel aa 


(vetween n and n+) j. 


Ko= mw ky K, = =->{h 


st 


Dl Lh Beale | 
G, a Fn G, ae G, | 


\K- = Torsional rigidity of wing in lb. inches per radian. 








AYALYSIS 


Because of the nature of this investigation and the difficulties 
involved in the neasurenent ae cia effect of one factor at a time 
on the torsional vibration of a wing, no experinental work was under- 
taken. Instead, the authors approached the problen from a ourely 
theoretical viewpoint, and the investigation was performed entirely 
on that basis. 

After a representative airplane for the investigation had been 
selected, it was necessary first to obtain the fol lewing infornation 
concerning it: 

A. The wing (Table I), considerine its aass and the masses of 
all bodies either attached to it or stored within it as 
being concentrated at a nunber of stations along its span: 

(1) the distance of each station fron the wing root in 
inches, 

(2) the nass polar manent of inertia tT, about the elastic 
axis of the wing of the mass considered to be concen- 
trated at each station in lb-inches seconds squared, 
and 

(S$) the rigidity I in lb-inches per radian, or its 
reciprocal, of the wing in torsion between each sta- 
tion. 

B. The fuselage (Table II), considering its mass and the nasses 
of all. bodies either attached to it or stored within it as 


being concentrated at a number of stations along its length: 








(1) the distance of each station from the fuselage nose 
in inches, 
(2) the bending monent of inertia about a ale oll axis 
verpendicular to the longitudinal axis of the fuselage 
&, in inches to the fourth power, and 
(3) the total mass m, considered as concentrated at each 
station in lbs. seconds squared per inch. 
After receist of the required information for the wing, it was possible 
to calculate the natural frequencies of the wing in torsion for as 
many modes of vibration as were desired, considering the wings as being 
built in to a stiff fuselage with an extresely high nonent of inertia 
conpared with that of each station along the wing. This calculation 
was actually carried out for two modes of vibration following Holzer's 
method as outlined on pages 228 and 229 of Ref. 1, an example of which 
has been aopended to this paper as Table III with an explanation in- 
cluded in the appendix. The results of this calculation have been 
tabulated in Table ‘IV and plotted on Fig. 1, and show that the natural 
frequencies for the first two nodes as determined by this calculation 
are 33.67 and 71.70 radians per second resvectively. 
This canpleted the first ~chase of the investigatior; and with the 
required infornation for the fuselage then at hand, it was possible 
to proceed with the second phase, nanely, the calculation of the natu- 
ral frequencies of the wing in torsion for as many modes of vibration 
as were desired. considering the torsional monent applied at the root 


of the wines by a flexing fuselage. The vroblems imnediately con- 


io 








fronting the authors in this phase of the investigation were those 
of deteruining (1) the torsional moment voroduced at the roct of a@ 
wing by a flexing fuselage and (2) the method of coupling this moment 
into the wing at its root. 

For the solution of the first of these problems a method devel- 
oped in Ref. 2 for the antisymetric bending of wings was applied 
to the flexing fuselage, considerin; the fuselage to be zade uv of 
two indenendent beems extending in opposite directions fron the lo- 
cation of the elastic axis at the root of the wing. This method has 
the advantage of yielding innediately the bending monent at any var- 
tioular station along a cantilever bean and the slope of the beam at 
that station as linear functions of the normal displacement of the 
beam. Consequently, the procedure followed was, first, to calculate 
the bending moments and the slopes, at the location of the elastic 
axis at the root of the wing, of both the portion of the fuselage 
aft of this location and the portion of the fuselage forward of this 
location resultinre fron a wnit downward displacement of the extrene 
end of both the tail and the nose. The bending moment at the elastic 
axis and the slove at that location of the after portion of the fuse- 
lage were designated as M, and «, respectively, and of the forward 
portion of the fuselage as M, and «x, respectively. 

Next, ames we fuselage is actually a continuous structure 
throughout its length, its slope on either side of the elastic axis 
must equal the slope on the other side of the elastic axis. This 


leads to the result that, since the initial displacenents of both the 
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tail and the nose were taken tc oe oositive downwards, in orler for 
the slove forward of the elastic axis to equal that aft of the elastic 
axis, the slepe forward of the elastic axis %, must be multiplied by 
the ratio (- xy /o<). This saze result would have been obtained hed 
the initial displaceaent of the nose been multiplied by this ratio 
(~ “v/s ) ; and since the bending moment developed is a linear func- 
tion of the disvlacement of the free end, the ben ling moment produced 
at the elastic axis b. the forward portion of the Pike lhe M, should 
also be multiplied by this sane ratio. “Ye then have that, for the 
continuous fuselage, the bending manente at the elastic axis due to 
the after and forward vortions of the “uselaze are given by the ex- 
pressions, My and ( Ky foxy) x My respectively. However, these two 
econponents oppose one anotner; consequently, in order to detemine the 
total bending monent M. fram the fuselage to be coupled into the root 
of the wing, one must be subtracted fron the other. If the direction 
of M, is taken to be the positive direction, it may readily be seen 
then that M.=M,- (- %/o.) Mp = My + (%/a4)M\. And since the bend- 
ing monent fron the fuselage at the elastic axis enters the wing as 
a torsional monent, M. is the torsional monent produced at the root 
of the wing, the enowmt of it entering each side of the wing being 7M, , 
assuming sy metrical twisting of the wing. 

for the nethod of coupling this amonrent into the wing at its root, 
one side of the wing was considered as a free body in torsion with this 
torsional moaert of 2M. avplied at its root. Assuming arbitrary unit 


angular deflections of the wing tiv, the Holser's calculations uaade 








during the first phase of this investigation yielded the following in- 
fornation for each frequency selected: 

(1, The torsional monent developed within the wing at the first 
station outboard of the root due to the rotational inertie 
forees within the wing ( S15, froa Table III,j, and 

(2) Tae angle of twist developed at the reot of the wing (4, 
fran Table iii). 

since in whi cadeuldeden the torsional monent developed at any 

particular station and the angle of twist at that station are given 

as linear functions of the arbitrary anguler deredugn or the wing 
tip, any desired anyle of twist at the wing root can be obtained by 
proverl, adjusting the arbitrary angular deflection of the wing tiv. 
Since the win; can be considered to be built-in to the fuselage, its 
angle of twist et the root should equal the slope of the fuselase at 
the wing's elastic axis, and in onder to obtain this angle of twist 

at the root 1¢ is necessary to multiply the original ar>vitrary angular 
deflection of the wing tin by the ‘ratio (ox, /3,) . daving aultiplied 
the original arbilrary angular deflection of the wing tip by this 
ratio, it is then necessary to multivly the torsional moment devel- 
oped within the wing at the first station éatboard of the wing root 

by this ratio also. Hence, this moment is then found to equal 

(x, /2,)ZL,wA, , and adding th.s to the torsional monent applied at 

the root of the wing b, the fuselage, a re-idual torsicnal monent or 
torgue on the wing of M=4Mo+ Ce) ZL04, is found. This residual 


torque is then the additional applied menent required in order to force 
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the wing to vibrate at the assuned frequency. This value of the re- 
Sidual torque is then plotted against the — frequency; an’ the 
process is then repeated for cther assumed frequencies, one point 
on the plot being obtained for each assuned _—n A complete 
example of the calculation by this method for one assumed frequency 
has been appended to this paper as Tables V(a/, V(b), Vi(a,, Vilv, 
and VIi,-with a brief explanation of then included in the appendix. 
The results of this calculation have been tabulated in Table VIII. 
After a sufficient number of points have been obtained, a curve 
may be drawn through then as has been done in Fig. 1. Again the 
points at which this curve crosses the frequenoy axis deteraine the 
natural frequencies of torsional vibration for the wing, for at these 
points the residual torque becones gero, and hence the additional 
applied moment required to force the wing to vibrate at that frequen- 
ey also becones zero. Fron Tig. 1 it ean be seen trat the natural 
frequencies for tne first two modes as determined by this calculation 
are 51.42 and 43.65 ee er beret respectively. (lee Table IX 
for tabulation of final results.,. These are reductions of 6.68% and 
39.1% respectively fron the frequencies of the first two modes found 
in the first phase of this investigation. Accordingly, it may be con- 
cluded that, whereas the consideration of e flexing fuselage has a 
snall but appreciable effect on the frequency of the first mode of 
torsional vibration of the wing, it has a very deciied effect in 
lowering the frequency of the second node. Probably, considering the 


trend of the curves of Fic. 1, this sane effect is carried on in pro- 
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gression to subsequent nodes of vibration; hence it is the studied 
opinion of the authors that this effect should be considered in the 
‘calculation of the natural — frequencies of the wing. 

The closing phase of this investigation was the determination of 
the fuselase deflection curves for each of the two modes of vibration 
gibi! nol above. This was accomplished with facility fron the cal- 
culations involved in the determination of the bending monents and 
slopes of the forwari and the after portions of the fuselage in the 
second phase of this investi;ation. The deflection at any station of 
the fuselage is designated as —_—" and columns so headed may be 
found in both Tables V~ (a, and Vi(b,. Again, the values of Yn 
given in Table VI(a,; aust be aultiplied s, the ratio (-a,/o,) in 
order to ¢ive then the correct magnitude with respect to those given 
in Table VI{b). YTuselage deflections are tabulated in Table X. A 
plot of the values of Yen calculated for the two natural frequencies 
found in the second nhase of this investigation was made and has been 
appended to this report as Fig. 2. <A perusal of this figure will in- 
dicate thet the deflection curves for the i for the two odes 
of wing torsional vibration arr very similar, there being no reflex 
curvatures along the fuselage length in either case. In the first 
mode the nose deflection is about one-seventh that of the tail whereas 
in the scecnd mode it is alnost twice that of the tail, fran which it 
ean be seen quite readily that for a given deflection of the nose the 
fuselage curvature will be auch greater for the first mole than it 


will for the second. 
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The relative angular defleotions of the wing at each station 
along its span can be determined very readily by referring to the col- 
unns headed @ in Table III for the calculation for a rigid fuselage 
and in Table VII for the calculation for a flexing fuselage. These 
values must be nultiplied by the ratio (x, /8,) for each frequency 
selected in the calculation for a flexing fuselage, as has already 
been done in the determination of the residual torque acting on the 
Wing, in order to determine the actuel nagnitudes corresponding to 
a unit downward deflection of the tail. This has been done and the 
results for the two modes tabulated in Table XI and plotted in Fig. 3. 

Fig. 6 is a schematic illustration of the two modes of vibra- 
tion, assuning a wit downward deflection of the extrene tail in each 


Case. 


Piha 








CCNCLUS LORS 


In the case of the airplene investigated herein, the consideration 
of a flexing fuselage has a small but appreciable effect on the fre- 
quency of the first mode of torsional vibration of the wing, but it 
has a very decided effect in lowering ib Pq unit of the second node. 

The deflection curwes for the fuselage for the first two nodes 
of wine torsional vibration are very siniler, there being no reflex 
curvatures along the fuselage length in either case. However, for a 
given deflection of the nose, the, fuselage curvature will be nuch 
greater for the first mode than it will be for the second. 

It must be understood that the above conelusions apply only to 
the pertiouler airplane which has been investigated herein. This 
paper is not subnitted with the intent to show that effects of sinilar 
magnitude can be expected for all airplanes, but simply that the effect 
should be investigated with the thought in mind that it might vrove 


appreciable, varticularly in the case of higher modes. 
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APPENDIX 
Se 2 A 


I. CALCULATION FOR THE RIGID ri FLAGE 

The method used is outlined on pages 228 and 229 of ne. 1, and 
is known as FHolver's method. The wing data (See Table I.) furnished 
for the airplane in question assumed the mass monents of inertia LI, 
of the wing to be concentrated at 7 scanwise stations, the first and | 
last stations being located at the tip and the root respectively. 

The notation used is the same as that for the fuselage and is denon- 
strated in rig. 4. 

A positive (clinbing, pitching angle at the tip ( n=! j of 
one radian was assumed, ( a@-i8 and for a given frequeney, the 
inertia torque Z 182, was calculated for station n-; . Thie 
inertia toraue multiplied by the torsional flexibility for panel 
length le gave the amount of twist or tie reduction in angle 2 be- 
tween stations n=; and n=2 . This angle of twist was then sub- 
tracted fron @ to give 2, , the angular deflection at station 
n=2 . Knowing <~, , the inertia torque at station n= 2 21,07, 
was then calculated. 

The remainder of the table was conpleted in like manner until 
the residual inertia torque = Tee, (at the wing root) was found. 
This value of inertia torque was tabulated in Table IV and plotted 
against wo infFig. 1. The residual inertia torque is the shaking 
monent which would be required at the wing root to cause the wing to 


vibrate torsimally at the assuned frequency w . > tne ome 
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of a rigid fuselage and assuming syunetric torsional vioration, = Tee 
must equal zero at e natural torsional frequency of the wing. Consequent- 
ly these natural frequencies can be found by plotting Zs TW, against 

t) to deternine points of intersection with the axis as was done 

in Fig. 1. This calculetion was carried out for the first two nodes, 

the final results avvearing in Table IX under "Rigid Fuselage”. 


n=7 J 
A sample calculation for 214 appears as Table III. 


II. CALCULATION FOR THE FLEXIBLE FUSELAGE 

(a) Outline of Method Used; 

The method used here for the fuselage is that derived in Pef. 2 
for the antisymmetric bending of airplene wings. A brief reswne of 
the method follows herewith. 

Using the notation demonstrated in Fig. 4, the shear at any ste- 


tion nN is given by 
Ain 
Se Z M,a) : (1) 
and the bending monent at any atation 1 is given by 


VM = Z dy, 6%) | i’ 


The slope at any station n+! is given by 


Co = ao, Tt a M.Ymy, (3) 


n+) 


~ae 





snd the deflection at any station n+\ ise given by 
Gantt ~ ga i? LX, Sade, a M, om, (4) 


where 
Vp,= change in slope fram n+! to n due to a mit force at n. 
Vm,= change in slope fron n+! to n due to a unit moment at n. 
dp = change in defleetion from n+! to n due to a wit force at n. 
Om,= cnange in deflection fron ne to hm due to a unit awaonent at n. 
The ‘nethod of calculation of these ~nmraaecters is outlined in sec- 

tion II (dy of this appendix. 

.cubstituting the expressions for S, and M, fraa equations (1) 


and (2) into equations (S, and («) 


Asn-i 


2 oe 4 
Say — eae Yen Z My: a DV", ya mM, We (x; Xn) (5) 
azn i aznel 
Gee = Ye - Lea, a aie Z yu. 7 oO du, a mM. yy (X= Son) (6) 


At the end of the fuselage, assune 


x= 

a? 
so that “= ¢- Ve, fo Se Ce (Sa) 
to = \-loa, — de m =| + LF, - d wm, -L¢ on Me L ¢ (6a) 


Continuing with equations (5) and (€) in like manner yields 
x. = 4D ai f (7) 


= £0 « 








ue= Oo — kh, @ (8) 


where coefficients fh , f, , qa, » and k, are independent of g , 
one couplete set being obtained for each frequency. The method of 
determination of these coefficients is outlined in section II (b) of 
this aovendix. For the present, these corfficients are assumed to 
be known for any particular frequency 

For antisymactric bending, the case where the fuselage is being 
shaken by a shaking moment Mcosuwt about the elastic axis of the 
wing, the deflection at the elastic axis is zero u= 0 and froa 
equation (3, ® - Se | With this value of dh , all of the deflec- 


ke 


tions y, may be found by means of equation (8), as can the bending ° 


. 


An 


monent at any station n , ZF mary, (KX: =m) 
A) f 
In the vartiecular caleulation with which this paper is concerned, 
the two quantities desired are the bending nonent coming in fron the 


tail or nose M, , and the slove of the fuselage at the elastic axis 


Xp 
n =e-1 

azb-1 acbo-l ; 
= z m- alg. (xz -X,) 7 Za mS ki (4, > X,) (10) 
Putting mag. = G; (12) 
and mak, = K, (12) 

s:b-\ 
also (x, - i (13) 








42 b-1 S*w-\ 
a4 


hen Mo =F (Ge 4) - ZKZ 4) (14 ) 


s 


But a Le ZL) - (44 ta--2,) 4G, (0, +2, +5- Lyn) +> 1G Uy 


. Qe fUCG 4ee) + ICs, 46,46,)+--L, pa Soe G,.,) 


heb- S24 
= 3 [ 2. va G,| (25) 
ixvb:! gsab-l Aeb- S>4 
and similarly ee \K.¢ ops = p Zz a a K. | 
co M= SUL ZG) -g 202 KI (16 


Referring to Table VI (»} 


Sta 


Ac is given by column (2, and 


| 


7G, is given by column (6), the second summations occur 


Sz} - 


ring in colwx1rs (3) and (7) which colwans give 


' 4zb-l S<A 
aye = a luke ze K. | (17) 
and a -= [ f. z e (18) 


Fron this it is seen that 


M=G&@-K¢ (19) 


and fron equation (7) a= AP - f (20, 
where f, and f, are given on line n=& under columns (4) and (8) 
respectively. 


The mnethod was repeated for the nose using the seme frequency 


- 22 = 












(Table VI(a,, to obtain M, and a, , ani the four values thus de- 
veraineda were tabulated at the top of Table VII. 

with an assuned positive (dowmward, deflection of one inch et the 
tail the total :.onent introduced at the elustic axis b, the fuselage 
is given oy | 


Mo= («,/ax,)M, + M, 










end the slope of the fuselage at the elastic axis is x, . 


(by Deteraination ef Coefficients 
with the original assuaptions at the end of the fuselage « = Pp 
end ysl , fron equations (7) and (8, 7 


Xx = A. aay ys 
Gin a qtr KD 





Re fz 0 q.< k, 


II 
O 


fron equations (fa) and (6a) 
Oe 
4.-= 9.- 4¢ 
k= 2, 





are ovtained £,= | 
Substituting equations (7, and (&, into equations (5) and (6), 
= ° 4a 


azn 


4 ; { ae - a ees and i 
Ree pitien= Rp -f, - Oy, Z ae K.9)(x;-X,) (21) 





4 
Z m.(qi- Keb) - Avy, Z 


Fast oe aa Te kp a a o = Dd, = milgi-kd) y Sdn, Z ml 


Nn art 


kK, (Kj -x,) 





» 








By equating tec. containing and those not ecnteaining ¢ on te 
two sides of eae of equations (21; and (22), the following equations 


are cvtaineds 


Cquating coefficients of Pf 


fey ~ on +O Me Sk + Wy, Z ae KO) (23) 
azn Lan! 
ee _ k, 2 Le wide ra m,k, a day va mM; k. (X.- X,) (24, 


kquatins constant terns: 


fa = He WY Ve eS moa OVe4 BZ Mm, Gi (er (25) 
Ath a A=N-1 s 
Fina = ¢n a Me —e — Wide Z MG ee du, ye mM, qe (x: - te (20 4 


Using substitutions and relaticns developed in (ll,, (12,, (13), (15) 


and (15): 


BR, time EK + pel Z Kad. (27) 


nei 


k eS iG ai a, - de, K; - a = iQ. z K, | (26) 
> ioe Ve, ya G 4 Vw, Z ee 2 om (29) 


ean te 4 Lifnn - de, ZG: ~ dn, Ze Z G, | (30) 


n'ont 


aks 


All the coefficients f. seule , » @nd q, can be found by pro- 


nN 
gressive calculation with the aid of e table (such as Tables Vi(e) and 


VI{bj, based on these equations. 


(c} Procedure Used in Filling out Tebles Vi(a, and VI(b/ 


st 
t 








Paraneters Ve , Ym, Gr, and dy were calculated (see section 
(a, of this appendix) and written in the spaces indicated. Like- 
wise the values for Hh were entered in the tables. 

Then, for a given frequency, values of m,” were calculated and 
entered in the appropriate colunn. 


Wext, in line n=) , the following values were entered in 


columns (1), (4) and (5) respectively: 
k =o f=) dl 


and in line n=z under column (1): \e ae 

The table was then worked across fron left to right starting with 
line n=: then vroceeding with line n=-z2 etc., each stev being in- 
dieated in the colimn heading. 

n comouting values to enter in colwnns (1,, (3), (5) and (7), one 
must remenber to use information appearing in the preceeding line. 

The remainder of the steps are self explanatory, the desired 


quantities of the calculation being M, end X, . 


(d) Calculation of Paraneters Ve, » Vm, » de, » and dm. 

fhe fuselaze data received (Table II., indicated bending monents 
of inertia equal to sero at each end of the fuselage, but in order 
to more nearly approximate the probable monent of inertia distribution 
in the regions fron the extreme ends to the next stations inboard, 
a trevesoidal distribution over these regions was assuned in both 
cases with the end ordinates approximately half the value of the next 


ordinates inboard. The assimed values were 


“I, (TAIL) = S00 in* I, (NOSE) = 3,500 in* 


~ 6 « 








A fuselage stiffness curve (suoh as Fig. 5, could be drawn for 
the given and assumed (end, values of bending monent of inertia, where 
the bending monent of inertia is plotted against fuselage distance 
as abscissa, such that the areas between succeedin;, stations would be 
trapezoids. If I, is taken as a convenient refe~ence value of the 
bending nonent of inertia for the fuselage as a whole, then at any 


point between stations n and n+} 


TL = iiGe ee 


“sing the monent area rethod: 
: : Ak& a wl Q b,, gy” \ Je b 
‘ns eer ET, b, Igelo* 2S), = Ett yee] 82, 


o 


de 
oF oe ain Eee Oo: Pa. HE) ‘Te [b,~ aby. CazP)| (32) 


Ww 


" eas 


= 
s, - febigs -inf 





v La by + ane (G2$%) — a,b.) (38) 


Referring to Tables V(a, and (b/, EL, was taken as (9'° . Colunns 
(1) and (£) were filled ir with values of Q. end gag, fron the data 


furnished. The value (a,+6,) in column (3) is determined fron 


a = @ +887 


n+\ 


Due to lack of availability of a seven nlace table of natural 
logarithns, conmucn logarithns were used and values converted to 


natural logerithus in columns (8) and (9) by the relation 


As. (ant On) = 2.302585 fog, (24) 


The remainijer of the table is self explanatory and follows froa 
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equations (51), ($2) and (35). The desired parameters appear in col- 


unns (8), (12) and (18). 


(e) Trocedure "sed in Filling out Table VII. 

‘gain using the sene frequency as was used in Tables Vita) and 
VIi(bj, Holeer's ecaleulation (es explained ir section I of this appen- - 
dix) was repeated in Table VII, the desired quantities being the 
angular deflection of the wing at the root , and the inertia torque 
at the next station outboard fron tne rcot = 1,88, - This inertia 
torque adjusted sc as to take the angular deflection of the wing at 
the root equal to the slope of the fuselare at the elastic — is 
(= 1,04} (u/6)  , and this adjusted ine-tie torque aided to half 
the honding moment introduced at the elastic axis vy the fuselace 
x M. gives the residual torque actin; at the wing root which would 
be required to nake the wine vibrate torsionally at the chosen frequency. 
MakM. + (Zz tow} (%/3,) 

At a natural frequency of the systen tris residual torque is zero. 
Table VIII gives a tabulation o! the results of this calculation. 

Fig. 1 shows a nlot of this residual torcue against © and the natural 
frequencies (first and second nodes) occur when this curve crosses the 
w axis, The final results are tabulated in-Table 1X under "i lexible 


kuselace”. 
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‘TABLF VIIT - RESULTS OF FLEXIBLE FUSELAGE CALCUL‘T ON 


CN 








RIGID FLEXIBLE 
FUSELAGE FUSELAGE PERCENTAGE 
(RADIANS /SEC) 


DECREASE 





TASLE IX ~- FINAL RESULTS 





' 
Crs 
er 

’ 
















SECONO MODE wW=43.6463057 


DEFLECTION 
FOR NOSE 
y'y (-0%5 foc) 
FOR TAIL 
4 a $ 
rr [re [ieoove | oiseits | Lovee | reness 
|| 7450] o7isii9 | otosii7 | 0.e646e4 | 1.240538 
0.030421 0.279369 
S| ° | 28856] 0.117336 | +O017152 | 0.078781 | +0\42743 
refelsaolo | 0 
ra[assol | o000n87 | —SCSC*~*~rtC mw 
438,00 - 0.017684 


+ 0.009384 
0.209373 
0.501356 


FIRST MODE ®=3).42148 


DEFLECTION 
FOR NOSE 
gx (©, /ox!) 


FOR TAIL 


DISTANCE 















STATION FROM 


n 























FUSE- 






LAGE 









NOSE 




































0.636222 






0.857362 
1.000000 





TABIO X = FUSELAGE DGPS CTIQg: 













FIRST MODE ®=31.42148 |SECOND MODE wW=43.64¢63057 


ANGULAR ANGULAR 
DEFLECTION 
(ADSUSTED) 
fF JB (%b //3_)X10" - 


/3 x (% /B,) x \O" 
| 6 | 133 | 03392151 | 0.388 |-0.247594 | -0.102 


NUMBER | DISTANCE 


















OF 





OF 
DEFLECTION ~- 


(ADI USTED) 







STATION 





WING 


FROM 





STATION 
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TABLE AI - AIG ANGULAR DEFLECTIONS 
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FIG. 5 - SAYPLT FUSELAGH STIFFNESS CUPVE AND NOTATICN USED 
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